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ABSTRACT 


An extremely accurate method for the measurement of mass and the 
lateral center of mass of a spacecraft has been developed. The method was 
needed for the Voyager spacecraft mission requirement which limited the uncer- 
tainty in the knowledge of lateral center of mass of the spacecraft system 
weighing 750 kg to be less than 1.0 mm (0.04 in.). 

The method consists of using three load cells symmetrically located at 
120° apart on a turntable with respect to the vertical axis of the spacecraft 
and making six measurements for each load cell. These six measurements are 
taken by cyclic rotations of the load cell turntable and of the spacecraft, about 
the vertical axis of the measurement fixture. This method eliminates all 
alignment, leveling, and load cell calibration errors for the lateral center 
of mass determination, and permits a statistical best fit of the measurement 
data. An associated data reduction computer program called MASCM has been 
written to implement this method and has been used for the Voyager spacecraft. 

The accuracy obtained was 0.14 mm (0.006 in.) for the center of mass and 0.01% 
for the mass determinations. 

The method and its associated data reduction computer program have been so 
designed that they are suitable for any other mass and center of mass 


determinations . 


1. SUMMARY 


An extremely accurate method for the measurement of mass and center of 
mass of a spacecraft is reported here. The method consists of using three load 
cells symmetrically located 120° apart on a turntable with respect to the 
vertical axis of the spacecraft and making six measurements for each load cell. 
These six measurements are taken by cyclic rotations of the turntable and one 
rotation of the spacecraft about its vertical axis* This method eliminates 
all alignment, leveling and load cells calibration errors for measurement of 
the center of mass, and permits a statistical best fit of the data, 

A data reduction computer program called MASCM has been written to imple- 
ment the method and has been used on the 750 kg Voyager spacecraft. The accuracy 
(standard deviation) obtained was 0.14 mm for the center of mass and 0,01% for 
the mass. 

2. METHOD 

The spacecraft is supported at the apices, 1, 2, and 3 of an equilateral 
triangle through an adapter ring and a load plate as shown in Figs. 1 and 2. 

Three load cells connect these three points to a common turntable which in turn 
rests on the ground. Two steps representing six sets of measurements are needed 
to eliminate alignment, spacecraft tiltings and load cell calibration errors. 
Details of the method are found in Appendix F. 

2.1 Center of Mass Measurement 

Step 1. 120 ° turntable rotations — — Three sets of measurements are taken 

with the load cell turntable assembly rotated 120° each time with respect to 
the spacecraft assembly and ground. Using measurement data for these three 
positions eliminates the load cell calibration errors (Appendix F) and the 
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Fig. 1 System Geometry and Notation 
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nonparallelism of top and bottom planes of the load cell turntable, see angle 


3 of Fig. 3 in the center of mass calculation. The load cell readings R. . are 

3-3 


different for each of the three positions, and generate nine readings tabulated 


as follows: 


^Measurement si 


I Load Cells 


where R.. “ load cell reading; i - load cell identification number and 
13 


j - load cell location relative to spacecraft points 1, 2, 3. 


For example, is the reading of load cell No. 3 at location 2 on the space 


craft. The coordinates x ^ , y. of the projection of the center of mass of the 


spacecraft assembly, including tare, on the interface plane P between the 


load plate and the load cell (Fig. 3) are calculated by (see Appendix F) : 














are the masses for each position, nominally different from each other because of 
calibration errors. Then, the coordinates X and Y of the projection of 
the center of mass are the mean value of the three sets of measurements. 

Ji = j <*i* + *2' + V 5 3115 ^ f ( y 1 + y 2 + y 3 ) C 5 ) 

Step 2* 180° Spacecraft rotation — Measurements are taken for the space- 

craft rotated 180° with respect to the load base plate held fixed relative to 
the ground. The turntable is then again rotated three times as before to 
obtain the second set of data. 

Rotating the spacecraft 180° eliminates the ground plane tilt error 
(angle a), load plate and fixture ring and concentricity errors (see Appendix C) , 

tt TT 

Using the same equations as above, a second set of coordinates X and Y is 
obtained. Then the coordinates and Y^^ of the center of mass of the space- 

craft with tare are: 


^OT = f < X ' + X ") and Y T0T = | < Y ’ + Y "> <6) 

2,2 Mass Measurement 

The mass of the spacecraft with tare W^q t Is the average value of the 
mass for each set 


w .1 y x 

TOT 6 A i 


R. 


iJ 


where are the load cell calibration factors. 


(7) 
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2.3 Tare 


The center of mass and mass of the tare must be removed from ^jOT 

and W^ot* Since the tare has a well-defined geometry, the coordinates 
and Y^ of the tare center of mass can be calculated with good accuracy. Hence 
a measurement as described above for the tare is not necessary but can be per- 
formed if desired. Also, since the tare is much lighter than the spacecraft, 

X T and Y t have only a small influence on the result. 

2.4 Net Spacecraft 

The center of mass coordinates X g ^ c and Y g ^ c and the mass of the space- 

craft are computed by 


x s/c 

w 

TOT 

' w s/c 

^TOT 

C8) 

Vc 

w 

w TOT 

w s/c 

Y 

TOT 

(9) 

w s/c 

= W TOT 

-w T 

(10) 


where W^, is the mass of the tare. 


2.5 Error Estimate 


The following can be used to estimate an upper bound of the errors before 
and after the measurement. 

a) Error on center of mass location of spacecraft plus tare 

The errors not eliminated by load cell and spacecraft rotations are 

estimated from the formula of Appendix F to which an origin of coordinates 

error Ax , Ay is added, 
o' o 


“ s4 | + | r ^ +&t o 


C11) 




( 12 ) 


— 
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Note that there is no error from the calibration factor because of the 
clocking of the load cells which measure the three reactions successively. 
The readout error AP is a random error made of three parts AP^, 
and APg (Appendix H) : 

AP^ is the readout drift error 

AP 2 is the readout error due to resolution 

APg is the readout fluctuation due to environment . 


b) Error on center of mass location of spacecraft 

From Eqs. (8) and (9) the center of mass error are 


AX 


W, 


s/c w, 


T0T AX 


S/C 


w r 


AY 


S/C W, 


TOT A? 


S/C 


(13) 


(14) 


where AX and AY are results of Eqs. (11) and (12) 


c) Error on mass of net spacecraft 

The error on the mass of the spacecraft is calculated by equation of AP 
of Appendix H 

AW = ^ W + 3 AP 


where the load cell calibration accuracy 


AX 


the 


being the error on the calibration 


is now introduced with 
factor A . 
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2.6 Statistical Errors 

The center o£ mass error estimates calculated in Section 2.5 are only 
bounds. After obtaining the measurement data, a best fit solution to the 

i f 

points , C 3c 2>5 r 2' > » in tlie ^ P lane ^ shown in Fig. 3 is 

obtained. The best fit criteria is that of parallel equal equilateral tri- 
angles and permit the determination of the standard deviation which is an 
indication of the quality of the measurements. Appendix G shows the derivation 
of the statistical best fit scheme and the resulting standard derivation. 

2. 7 Load Cell Calibration and Requirements 

Accurate load cell calibrations, as discussed in Appendix F, are needed 
only to obtain an accurate measurement of spacecraft mass but are not required 
for the center of mass. The ranges of needed accuracy for the Voyager space- 
craft measurement are from 45 to 90 kg and from 275 to 365 kg. Sample 
calibration curves are shown in Appendix D. 

The following were the instrumentation requirements to meet the Voyager 
spacecraft measurement accuracy limits: 


1 . 

Load cell range 

1000 kg 

2. 

Readout resolution 

.05 kg 

3. 

Load cell readout system 
calibration accuracy 

.1% 


f in 45-90 kg and \ 
\275-365 kg range/ 


4. 

Load cell location accuracy 

.05 mm 
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Data Reduction Computer Program MASCM 


3.1 Program Summary 

The computer program MASCM provides the means of computations of the 
mass and the center of mass of the spacecraft from data obtained by the 
measurement operations as described in Section 2, and the estimation of the 
standard derivations defined by Eqs. (2) through (9) of Appendix G. The 
usage of MASCM is described below. Relevant mathematical derivations and a 
listing of the program are in Section 2 and the Appendices. 

The lateral center of mass location (x, y) output from MASCM is referred to the 
spacecraft coordinate system rather than to the load cell coordinate system 
used in Section 2. The coordinate transformation from the load cell to 
the spacecraft coordinate systems is described in Appendix B. In addition 
Appendix B shows the manner in which the program handles the cases of tare 
measurements and spacecraft plus tare measurements, the method of correction 
for these measurements and the calibration data. It is also shown that the 
lateral center of mass of the tare is not required for the calculation of 
the spacecraft lateral center of mass. 

The measured mass and center of mass values of the Voyager spacecraft 
have been determined through the use of the MASCM program and are presented 
in Table 1. All uncertainties are la (one sigma) value. 


■-SS5S5 
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TABLE 1. Mass and Center of Mass Prediction For 

The Voyager Spacecraft (Based on Measurement Data) 


SPACECRAFT 

(partial) 

MASS 

kg 

(lb) 

X 

cm 

(in) 

Y 

cm 

(in) 

ACG 

cm 

(in) 

(PTM) 

706.266 ± 0.083 

+1.191 ± 0.010 

+2.332 ± 0.010 

+0.014 

VGR77-2 

(1557.05 ± 0.183) 

(+0.469 ± 0.004) 

(+0.918 ± 0.004) 

(±0.006) 

(S/C 31) 

687.814 ± 0.091 

-0.079 ± 0.034 

+1.143 ± 0.034 

±0.048 

VGR77-2 

(1516.37 ± 0.200) 

(-0.031 ± 0,013) 

(+0.450 ± 0.013) 

(±0.019) 

(S/C 32) 

709.350 ± 0.096 

+0.132 ± 0.034 

+3.625 ± 0.034 

±0.048 

VGR77-3 ' 

(1563.85 ± 0.212) 

(+0.052 ± 0.013) 

(+1.427 ± 0.013) 

(±0.019) 


3.2 Program Description 

The MAS CM program is a computer program written in the MBASIC language. 

It is designed to compute the mass and the center of mass and their associated 
standard deviations from sets of data using Eqs. (2) through (9) of Appendix G 
A listing of the program is shown in Appendix A. 

The MASCM program is called either from the MBASIC file or from tape 
and is loaded into the computer for execution. It can be executed either from 
the demand terminal or through the batch processor. For a batch run the input 
data are punched on cards and the deck of cards submitted to the computer. Fo 
the demand terminal, the input data are entered stepwise into the computer as 























the computer requests them (see Appendix E, Sample Problem). The order in 
which the data are punched on cards or entered on the demand terminal is the 
same for both batch and demand terminal modes. 

The technique as described in Section 2 requires measurements to be 
performed with the rotation of the load cell turntable assembly to three 
different positions and with the spacecraft in each of the two positions, 

0° and 180° rotation with respect to the tare. Therefore, two tables of load 
cell data for the spacecraft plus tare case representing 18 inputs are needed 
for inputs. 

Because of design simplification of the spacecraft handling scheme during 
measurement, there is only one set of load cell accommodation holes on the 
load plate (see Fig. 2 of Section 2) . Therefore, the load plate does not have 
the capability of being measured in the two positions and only one table of 
load cell data for the tare is obtained. The MASCM program is designed to 
have the capability of handling cases of either one table (as the tare case) or 
two tables (as the spacecraft plus tare case) of load cell data. The particular 
case desired for computation is implemented by exercising the option — Type 1 
or Type 2 respectively — in the program. 

The input can be either in the metric units or the U.S. customary units. 

3.3 Program Usage 

3.3,1 The input data to the program consists of: 

. The test (measurement) identification number 
. The radius, r of the load cell location, see Fig. 1 and 
Eqs. (3) through (6) of Appendix B 
. The angle, fl of the load cell system orientation, see Fig. 1 
and Eqs. (7) and (8) of Appendix B 

original page IS ii 
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. The coordinates x q , y Q of load cell system origin 0^, 
see Fig. 1 and Eqs. CH) and (12) of Appendix B. 

. The load cell readout data. R, . as tabulated in Section 2 

ij 

. The configuration of the spacecraft that is undergoing 
measurement. 

. The tare mass and CM data 

. Option (Type 1 or Type 2) to select the case of one or two tables 
of load cell readings 

. Option (1 or 2) to select input units (metric or U. S. customary). 


3.3.2 The key to the Input/Output Parameter: 


RADIUS: Average radius of load cell locations 

MASS: Mass of a specified portion of spacecraft 


XBAR: Abscissa of a CM location 

YBAR: Ordinate of a CM location 

SIGMA: Standard deviation of a CM location 

SGMAW: Standard deviation of mass 

SGMCAL: Standard deviation from the load cell calibration curves 

R1(I, J) 9 R2(I,J) refer to the two sets of load cell readings 
corresponding to the 0° and 180° rotations of the spacecraft. 

Each set is as shown in Section 2. These values must be entered 
by row. An alphabet T following any of the above parameters 
indicates that the tare is being referred to, for example, 

XBART = abscissa of the center of mass location of the tare alone 


3.3.3. Sample Runs tr earns 


The user may store this program on a magnetic tape. For example, 

for tape Y000, the runstream is as follows: 

@RUN X£X,XXXXXX,XXXXX 
@CAT,P FILE. 

@ASG,A FILE. 

@ASG,T T,T, Y000R 
^REWIND T. 

@C0PY,G T. ,FILE. 

@FREE T. 

QMBASIC (Accesses MBASIC) 

LOAD T FILE.MA.SCM' (Accesses MASCM program for execution) 

RUN (Execution of MAS CM program) 

\ 


I 

(-INPUT DATA (demand terminal entries or batch input deck) 


Output print0ut 

I 

EXIT 

0DELETE FILE. 

0FIN 


The user may store this program in the MBASIC file under his pass- 
word if he wishes to access the program from the MBASIC file for 
subsequent runs. In this case the runstream is as follows: 


■DfeftP, IS 
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@RUN XXK t mxKX,SXXXX (Run card) 

0MBASIC XXXXXX (Accesses MBASIC and its files) 

LOAD 'MASCM* (Accesses MASCM program for execution) 

RUN (Execution of MASCM program) 


— Input Data (demand terminal entries or batch Input deck) 


| — Output Printout 


The order of the input data is shown in the Appendix E (Sample 
Problem) . This order is valid for both demand terminal and batch 
runs as discussed in the beginning of Section 3,2* 

4. COORDINATE TRANSFORMATION 

Define x , y as the apparent center of mass of the spacecraft being measured 
P P 

in the load cell system of coordinates, then a corresponding center of mass x, y 
in the spacecraft system is given by Eq. (10a) of Appendix B and repeated here 


for convenience. 


cos0 sin0| f x 


-sin0 cos0 j^pj 


where x Q ,y o is the origin of load cell system with respect to spacecraft 
system and Q is angle of rotation of the spacecraft system with respect to 
the load cell system (Fig. 4 ) . 
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5. LOAD CELL CALIBRATION AND READOUT CORRECTION 


All load cell data R. . obtained during the mass and center of mass 
measurement of the spacecraft must be recorded in accordance with the require- 
ment delineated in Section 2.1. The readouts R. , must be corrected for the 

ij 

load cell calibration error before input to the program. The calibration 
error is obtained in Appendix D. An average value for each corrected R„ 
must be generated if more than one complete set of nine measurements is 
performed since each complete set must contain only nine . 

The value of SGMCAL (a^^) used for standard deviation estimate in the 
program is the average value of the standard deviations of all three load cells 
throughout the applicable load range. It Is established from data obtained by 
the load cell calibration test. Typical calibration data are presented in 
Appendix D. 


6. REMOVAL OF TARE 

As stated in Section 2 the load plate remains in the original position 
after the first set of measurements while the spacecraft is rotated 180° for the 
second set of measurements. To obtain the net spacecraft center of mass in 
the spacecraft coordinate system, the following mathematical expressions have 
been derived and are included in the MASCM program whe-re is the mass of 

the load plate plus the mass of the adapter ring (tare). The derivation 
Is shown in Appendix G. 


- / w s/c + w t\ - - 

3C . = I l I X . y 


s/c \ w, 


s/c 


w s/c + ^ 


s/c 


w, 


s/c 


' w s/c + ^ 


s/c \ w 


s/c 


V ^ 


2 + rr 2 

W S/C+T W T 
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Appendix A. PROGRAM LISTING 
MAS CM 


£ print *input test no.’ 

3 INPUT TESTNQ 

5 PRINT ’INPUT 1 IF UNITS ARE METRIC » 2 IF UNITS ARE U.S. CUSTOMARY’ 

7 INPUT SYSTEM 

9 REAL X ( 6) » Y (6) pSUMX£7) p SUMY (7) t A (3) ’B (3) pDC4) rDl (3 ) pD2(3) pD3£3> »D4<3) 

11 REAL R(3p3>?R1£3p3)pR2{3p3) pW(3)pW1<3) pW2£3)pXD(6)pYd(6) 

13 PRINT ’INPUT THE CONFIGURATION* 

15 INPUT CONFG$ 

17 PRINT ’INPUT RADIUS? WEIGHTTpXBART » YBART rSGMAWT »SIGMAT»SGMCAL» 

21 INPUT SR ? WT ? XBART ? YBART ? SGMAWT ? SGMAT ? SGMCAL 
23 PRINT ’INPUT XO?YO? THETA IN DG.* 

25 INPUT XOrYO? THETA 

27 PRINT ’INPUT R1(I?J) ROW BY ROW’ 

29 INPUT Rl(lrl) ?R1£1?2) ?Rl(l?3) ?R1£2?2) ?Rl(2?3) » R1 <2?1) ? R1 C3 p 3) ?R1£3?1)»& 
Rl(3?2) 

30 PRINT’ INPUT 1 OR 2 FOR ONE OR TWO TABLES OF L/C READINGS’ 

31 INPUT MESUD 

32 XBART2=X8ART 

33 YBART2=YBART 

34 SGMAT2=SGMAT 

35 IF MESUD=1 THEN GO TO 42 

36 XBART=0. 

38 YbART-G . 

40 SGMAT=0. 

42 R=Rl 

44 GOSUB 500 
46 W1=W 

48 COSTH=COSC THETA) 

50 SINTH=SIN( THETA) 

51 XCl)=XO+COSTH*XD(I)+SINTH*YDf I) FOR 1=1 TO 3 

52 YU)=YO-SINTH*XO(I>+COSTH*YD( I) FOR 1=1 TO 3 

53 IF MESU[)=1 THEN GO TO 62 

54 PRINT ’INPUT R2£I?J) ROW BY ROW* 

55 INPUT R2<lrl) »R2£1 p2) pR 2(1?3) ?R2(2»2) ?R2(2r3) ?R2(2?1) »R2(3?3) 
R2f3?l)?R2£3?2) 

56 R=R2 

57 GOSUB 500 

60 X(I+3)=-XO-COSTH*XDtI)-5INTH*YD(I) FOR 1=1 TO 3 

61 Y(l+3)=-YO+SINTH*XD<I)-COSTH*YD(I) FOR 1=1 TO 3 

62 W2=«l 

63 IF MESUD=2 THEN GO TO 67 

65 X(I+3)=X(I) FOR 1=1 TO 3 

66 Y(I+3)=Y(I) FOR 1=1 TO 3 

67 SUi-1X C 1 ) =0 ? SUMY ( 1 ) =0 
70 FOR r=l TO 6 

73 SUMX 1 1+1 ) =SUMX ( I ) +X ( I ) 

75 SUMY C 1+1 ) =SUMY { I ) +Y ( I ) 

80 NEXT I 

85 XBAR=SUMX(7)/6 

90 YBAR=SUMY(7)/6 

95 ALFA= (2*SUNX (4)-SUMX (7) )/6 

100 BETA=C2*SUMY(4)-SUMY(7) )/6 

110 BB=(3*YC1)+3*YC4)-SUMY(7)-SQRC3)*CXC2)-X£3)+X(5)-X(6)) )/12 
120 AA=(3*X(1)+3*X(4)-SUMX(7)+S0R{3)*(Y(2)-YC3)+Y(5)-Y(6) ) )/12 
130 A £ 1 ) =AA 

132 A ( 2 ) =- £ A A+SQR £ 3 ) *BB ) /2 
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134 A(3)=-(AA-SQR(3)*BB)/2 
136 B ( 1)=BB 

138 B(2)=(SQR(3)*AA“BB)/2 
140 B(3)=-(SQR(3)*AA+BB)/2 
143 0(1) =0* WEIGHT=0 
145 FOR 1=1 TO 3 

147 WElGHT=WEl6HT+( W1 ( I ) +W2 ( I ) )/6 
150 01 ( I ) =X ( I ) - ( XBAR+ALFA+A ( I ) ) 

162 D2(I)=Y(Y)-(YBAR+BETA+B(I) ) 

164 D3 ( I ) =X ( 1+3 ) - ( XBAR-ALFA+A ( I ) ) 

166 D4 ( I ) =Y ( 1+3 ) - ( YBAR-BETA+8 (I) ) 

190 D ( 1+1 ) =D1 ( I ) **2+D2 ( I ) **2+03 < I ) **2+04 ( I ) **2+D ( I ) 

200 NEXT I 

210 SlGMA=SQR(D(4)/5) 

212 5GMWSQ=0 

213 FOR 1=1 TO 3 

215 SGMWSQ= (WEIGHT-Wl ( I ) ) **2+ ( WEIGHT-W2 ( I ) )**2+SGMWSQ 
217 NEXT I 

220 SGMAW=SQR (SGMWSQ/5+SGMCAL**2 ) 

240 WSC=WEIGHT-WT 

250 X3ARSC=WEIGHT*XBAR/WSC-WT*XRART/WSC 
260 YBARSC=WEIGHT*YBAR/WSC-WT*YRART/WSC 
270 SGMASC=SQR( (WEIGHT*SIGMA/WSC)**2+(WT*SGMAT/WSC)**2) 

275 S6MWSC=SQR ( SGMAW**2+SGMAWT**2 ) 

276 PRINTSSS 

277 PRINT USING’ TEST NO. SS*9S» STESTNO 


MEASURED 


WEIGHT' 


MEASURED 


MASS* 


CONFIGURATION: 


:confg$ 


INPUT* 


THEN GO TO 295 

WEIGHT READINGS (LB>» 


READINGS (KG)’ 


277 PRINT USING’ TEST NO. SS**»:TESTN0 

278 PRINT\ 

280 IF SYSTEM=1 THEN GO TO 284 

281 PRINT * 1, MEASURED C.G. AND WEIGHT* 

282 GO TO 285 

284 PRINT » 1. MEASURED C.M. AND MASS* 

285 PRINTS 

286 PRINT USING • CONFIGURATION: «':C0NFg$ 

£87 PRINTS 

288 PRINT* INPUT* 

289 PRINT 

290 IF SYSTEM=1 THEN GO TO 295 

292 PRINT * WEIGHT READINGS (LB) ’ 

294 GO TO 298 

295 PRINT * MASS READINGS (KG)’ 

298 FMT5$= * ( /3 ( 3B35K . 355 >) ’ 

300 FMT1S=* (/2(3(3B355.3SS)6B) ) * 

302 IF MESUD=2 THEN GO TO 310 

305 PRINT USING FMT5S : R1 (1 , 1 ) »R1 ( 1 »2> rRl { 1 , 3J ,& 

Rl(2r2) »R1C2»3) pRl(2rl)»Rl(3p3) f Rl(3»l)pRl(3p2) 

308 IF MESUD=1 THEN GO TO 313 

310 PRINTS* CASE 1. CASE 2.’ 

312 PRINT USING FMT1$:& 

Rl(l»l)»Rl(lr2) rRl(l»3) »R2(1»1) »R2(lr2) »R2(1*3) >8 
Rl(2»2)»Rl(2f3),Rl(2»l) »R2(2»2 ) pR2(2*3) pR2(2»1)pS 
Rl(3r3) * Rl (3 r 1 ) »Rl(3»2I»R2(3f 3 ) p R 2 ( 3r 1 ) r R2 ( 3? 2) 

313 FMT2$=»<(9B*XBAR (IN) : *556.4%/9B’YBAR (IN) i *5X.4«/8B* SIGMA (IN)!*8 
5SS . 4SS/ ) > ’ 

314 FMT22$=’<(9B*XBAR (CM) : *5SS.495/9B*YBAR (CM) : * 5SS, 4S8/8R ’SIGMA (CM):’& 
55U4%/)>* 

316 FMT7$=*<(9B»XBAR (CM) : *5%.4»16B’XBAR ( IN) : ’ 5S6.4S5/98 ’ YBAR (CM): '8 
5&.4*16B'YBAR ( IN) : ’ 535.4*) >’ 


CASE 2.’ 


(IN) : 


317 FMT77$='<{/aB* SIGMA (CM) i » 5%. 4S615B* SIGMA ( IN) I * 5SS.4«> >» 

318 FMT3S='<(7B*WEIGHT (LB) : »5%.4K/8B' SGMAW (LB) ! * 5%. 4*/) >« 

319 FmT33S='<C7B» MASS (KG) J »5«.4K/8B'5GMAW (KG) : * 5*. 4*/) >♦ 

320 FMT8$=*<(9B*XBAR (IN) : *5K.4(K16B*XBAR (CM) : *5?S.4fc/9B* YBAR (IN)J*a 
5$U4K16B'YBAR (CM) I »55B.4K/)>* 

321 FMT80$=»<(/8B* SIGMA ( IN) ! »5S.4*15B»5IGMA (CM) ! »55K.495)>« 

324 FMT4$=K(7B»SGMCAL (LB) J *585. 4*/) >» 

325 FMT6$=*<(7B*SGMCAL (KG) ! * 5S .4*1/) >' 

326 PRINT\ 

327 IF SYSTEM=1 THEN GO TO 332 

328 PRINT USING* RADIUS (IN)! 85%. 855885%* I SR 

329 PRINT USING FMT4SSSGMCAL 

330 GO TO 336 

332 PRINT USING’ RADIUS (CM): SfiSS.***** JSR 

334 PRINT USING FMT6$:SGMCAL 
336 PRINT 
338 GOSUB 600 

340 PRINT* OUTPUT* 

342 PRINT 

344 IF SYSTEM-1 THEN GO TO 353 

345 PRINT USING FMT 8$ ‘ X8AR t XBAR2 * YBAR r YB AR2 

346 PRINT USING FMT88S:SIGMA»SIGMA2 

347 FMT9$=*<(/7B* WEIGHT (LB) ! »5».4*15B* MASS (KG) ! ’ 5SU4Sfi)>* 

348 FMT99S=*<(/8B*SGMAW (LB) i » 5%. 49515B * SGmAW (KG) ! *5%.485/)>* 

350 PRINT USING FMTGSJWEIGHT* WT2 

351 PRINT USING FMT99SJ SGMAW »MAW2 

352 GO TO 359 

353 PRINT USING FMT7$ S X8AR t XBAR2 » YBAR r Yf3 AR2 

354 PRINT USING FMT77S! SIGMA »SIGMA2 

355 FMT10$=»<(/8B* MASS (KG) : »585.48514B* WEIGHT (LB) ! »5SG,4%)>* 

356 FMTli$=*<(/8B*SGMAW (KG) :« 5SS.4SS15B* SGMAW (LB ) ! * 585.4%) >* 

357 PRINT USING FMT1G$;WEXGHT* WT2 

358 PRINT USING FMT11$ISGMAW»MAW2 

359 PRINTW 

360 PRINT* 2, TARE* 

362 PRINT\ 

370 PRINT * INPUT* 

372 PRINT 

376 IF SYSTEM=1 THEN GO TO 386 

380 PRINT USING FMT2$!XBART2f YBART2fSGMAT2 

382 PRINT USING FMT3$!WT»SGMAWT 

384 GO TO 390 

386 PRINT USING FMT22$lXBART2>YBART2iSGMAT2 
388 PRINT USING FMT33$:WT»SGMAWT 
390 PRINTW 

395 IF SYSTEM=1 THEN GO TO 405 

400 PRINT USING* 3, MEASURED ft C.G, AND WEIGHT* SCONFG'S 

402 GO TO 410 

405 PRINT USING* 3. MEASURED 4 C.M. AND MASS*:CONFG$ 

410 PRINT 

415 PRINT * OUTPUT* 

420 PRINT 

425 IF SYSTEMS! THEN GO TO 438 

430 PRINT USING FMT8$!XBARSCrXSC2tYBARSC,YSC2 

431 PRINT USING FMT88$:SGMASCrSGSC2 

434 PRINT USING FMT99!WSCrWSC2 

435 PRINT USING FMT99S:SGMWSC»SMWSC2 
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436 GO TO 480 

438 PRINT USING FMT7$ I XB AR SC » X5C2» YB ARSC » YSC2 

439 PRINT USING FMT773>:SGMASCr SGSC2 

440 PRINT USING FMTIOS ; WSC > WSC2 

442 PRINT USING FMT11$*.SGMPSC rSMWSC2 
480 END 

•500 W(1)=R(1»1)+R(2»2)+R<3»3) r W{2)=R<lf £>+R(2»3>+RC3 f 1) 

510 W 1 3 ) — R ( 1 f 3 ) +R ( 2 f 1 ) +R 1 3 f 2 ) 

520 V=3*SR/2r U=SQRC3)*SR/2 

530 XuCl)=V*R(ia)/WCl) * XD(2)=V*RC3>1)/VJ(2) c XD {3)=V*R (2f 1) Al (3) 

540 YDCD=U*(R(2»2)"R(3»3) >/Wtl) » YD(2)=U*(R< 1*2)-R(2»3) )/W(2) 

550 Yj ( 3 1 = U* ( R ( 3 ► 2 ) “R ( 1 » 3 > ) /W ( 3 ) 

575 RETURN 

600 IF SYSTEM=2 THEN GO TO 630 

602 XbAR2=XBAR*.3937 

604 YBAR2=YBAR*.3937 

606 SXGMA2=SI6MA*«3937 

608 XSC2=XBARSC*.3937 

610 YSC2=YBAPSC*,3937 

612 SGSC2=SGNASC*.3937 

614 WT2=WEIGHT*2. 204622 

616 MAd2=SGMAW*2. 204622 

618 WSC2=WSC*2, 204622 

620 SNWSC2=SGMWSC*2.204622 

625 Go TO 650 

630 X8AR2=XBAR/,3937 

632 YtJAR2=YBAR/,3937 

634 SIGMA2=SIGMA/.3937 

636 XSC2=XBARSC/.3937 

638 YSC2=YBARSC/,3937 

640 SGSC2=SGMASC/.3937 

642 WT2=WEIGHT/2. 204622 

644 MAW2=SGMAW/2. 204622 

646 WSC2=WSC/2. 204622 

648 SMWSC2=SGMWSC/2.204622 

650 RETURN 
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Appendix B. COORDINATE TRANSFORMATION 


Figure B.l shows the position or location of the load cells on the 
fixture load plate (tare). Points (l), (_2) , and ^3) are the load cell locations 
which, nominally, form an equilateral triangle. The coordinates of the vertices 
referred to spacecraft system are indicated on the figure. For the convenience 
of describing the coordinate transformation, numerical expressions corresponding 
to Voyager have been included. Point G is the centroid of the triangle and 
its coordinates (x_, y„) are: 

x^ “ (x^ + x 2 + /3 = 0.00630 in. Cl) 


y G “ ^ y l + y 2 + y 3^ 3 = °* 00267 in * 


( 2 ) 


where (x^) , (x^) and(x 3 y 3 ) are the coordinates of points (l),(2) and (|) 
respectively. 

The calculated radii between G and points (T) , ( 2 ) and (3) are: 


r i 0 1 < x i - v 2 + (y i 

r 2 = j^ x 2 ” X G^ + ^ y 2 

r 3 " [C x 3 " X G )2 + (y 3 


- y G ) 2 ! 

1/2 

= 9.49879 

in. 

(3) 

“ y G )2 

1/2 

= 9.49908 

in. 

(4) 

“ y c )2 l 

1/2 

» 9.49917 

in* 

(5) 


These three radii would be identical and have a design value of 9.500 in. if 
triangle (T) ( 2 ) ( 3 ) was a perfect equilateral triangle. However, due to the 



Fig* B.l. Actual Dimensions of the Voyager Mass and CM 
Measurement Fixture Load Plate 
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manufacturing tolerance, it is only approximately so and the radius r is 
taken as the average of the three radii 


r + r„ + r 

r = _± 1 i = 9.49901 in. « 

Note that triangle (l) ( 2 ) (3) is equilateral within .001 in. 

Gall x , y the load cell coordinate system with its origin at 0 and 
P P P 

0 the angle of rotation of x, y system with respect to x^ , y^ system, then: 

. „ 7 1 ~ y G 4.791 - 0.00267 „ . 

sxn 9 = 9 1 4990 1 -0.499673 xn. (, 

„ *1 " X G 8.2342 - 0.00630 n . ,, 

cos 0 = ,, .... = 0.866214 in. (.1 

r 9.49901 


Eqs. (7) and (8) imply that 0 = -29.9784 . 

Fig. B.2 depicts the coordinate transformation for any point OCx^, y^) in 
Xp* y p system to Q(x,y) in the x,y system. 

*.* x=0F=0D+DE+EF 


x. e. , 


y = OC = 0A + AB + BC 


x = x + x cos 0 + y sin 0 
o p J p 


y = y + y cos 0 - x sin 0 
3 o J P p 
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where x , y are the coordinates of 0 relative to x,y system, 
o J o p 

Rewriting Eqs. (9) and (10) in matrix form, we have 



Fig. B.2. Coordinate Transformation 

Eq. (10a) is now written for the Voyager geometry. The coordinates of 0 
referred to a system parallel to x^, y^ with origin at G is (-GO^, 0). By 
using Eqs. (9) and (10), the coordinates of point 0^ with respect to the x, y 
system are given as: 
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( 11 ) 


x = (-GO ) cos 0 + x„ = -4.107788 in. 
o p G 

y = (-GO ) sin 0 + y„ = -2.370499 in. 
o p « 


( 12 ) 


— 1 

where GO = •=■ r 
P 2 

i i 

Substituting the numerical values of x^, and 9 and calling x , y 

the coordinate corresponding to the 0° position of the spacecraft we obtain 
the transformation 


x ? = -4.107788 + 0.866214 (x ) -0.499673 (y ) (13) 

P P 

y* = -2.370499 + 0.499673 (x p ) +0.866214 (y p ) (14) 

For the second se 1 " of measurements the spacecraft with the adapter ring 
was rotated 180° with respect to the load plate or tare (See Sec. 2). Recalling 
that the spacecraft coordinate system is fixed with the spacecraft, the 
appropriate transformation, corresponding to the 180° position of the space- 
craft, is obtained from Eqs. (13) and (14) by a change of sign. 

x" = 4.107788 - 0.866214 (x p ) + 0.499673 (y ) (15) 

y" = 2.370499 - 0,499673 (x ) - 0.866214 (y ) (16) 

P P 

Eqs. (13) through (16) are the transformation equations based on load 
cell geometry as shown in Fig. B.l. 
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Appendix C. REMOVAL OF THE TARE 


From Fig. C-l the spacecraft is at 0° position relative to the tare, 

then, 

W T x, r + W g/C (a + x g/c ) cos a = (W s / C + W T ) x cos a (1) 

where 

t 

x = combined spacecraft and tare C.G. 

a = nonparallelism of the top of load plate and the referenced 
horizontal ground plane. 

From Fig. C-2 the spacecraft plus the attached adapter ring only is at 
180° position relative to the tare, then, 


“Vt " W s/C (a “ x s/C 3 COS “ = (W S/C + V X COS ° 

where 

TT 

x = combined spacecraft and tare C*G. 

From Eq, (1) and (2) above 


( 2 ) 


S/C S/C 


Similarly 


cos a 

- (W S/C + V ' 



w s/c + 

W T 

(i 

x s/c 

H s/c 


\ 



_ H s/c + 

W T 



x s/c 

w s/c 


X 


w s/c + 

W T 

m 

y s/c 

w s/c 


y 


it 


(3) 

(4) 


where x, and y are the C.G. of spacecraft plus tare and are not a function 
of a. 




IS 
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Eqs. (3) and (4) show that the C.G. coordinates of the tare are not 
required in the calculation of the spacecraft's C.G. Only the weight of 
the tare W,p is a required parameter. Then the is in the form of 


a s/c 



s/c + w t\ 
w s/c / 


cr 


( 5 ) 


where a is calculated standard deviation on the spacecraft plus tare C.G. 
coordinates based on input data. However, the standard deviation on the space- 
craft weight estimate remain the same as that shown in Eq. (9) of Appendix G. 



W, 


+ a* 


S/C+T 


w. 


Tj 


1/2 


C6) 


where a TT and a T . are calculated standard deviation on the weight of the 

W S/C« W T 

spacecraft plus tare, and of the tare, respectively. 

As depicted in Fig. 2 of Section 2, the tare for Voyager comprised the 
adapter ring plus the load plate with a total weight of about 140 kg (310 lbs). 
The field joint is located at the interface of these two parts. The space- 
craft is detached from the load plate at the field joint for the 180° rotation 
with the adapter ring attached to it. Since first, the C.G. of the adapter 
ring is very close to the vertical axis of the spacecraft within 0.0025 cm 
(0.001 in.) and second, the weight of the ring is very small as compared to 
that of the spacecraft (1%) , the ring can be included in the load plate tare 
weight in Eqs. (3) and (4) . 


Appendix D. LOAD CELL CALIBRATION CURVES 



iiDaaao: 



Fig. D.l. Calibration Curve, Load Cell No. 1 (S/N 460), Yellow 
(Calibrated 1/30/76 and 2/2/76) 




CORRECTION FACTOR A, LB 


ADD "CORRECTION FACTOR" TO 
0.20- INDICATED READING TO OBTAIN 
TRUE READING 


a_ A , — 0.000 LB 
CAL 



- 0.10 


a cM _ = 0. 1668 LB 


0-200 LB RANGE 
A= 0 

450-900 LB RANGE 
A = -HJ.00082P -0.7335 


150 200 450 500 


550 600 650 700 750 800 850 900 

P LB 


Fig. D.2. Calibration. Curve, Load Cell Wo. 2 (S/N 461), Blue 

(Calibrated 3/17/76) 





LO 

ro 




P LB 


Fig. D.3. Calibration Curve, Load Cell Wo. 3 (S/N 462), Red 
(Calibrated 2/3/76 and 2/4/76) 



Appendix E. SAMPLE PROBLEMS 


INPUT TEST ND. 

71 

INPUT 1 IF UNITS RRE METRICS 2 IF UNITS RRE U.S. CUSTDMRRY 
?1 

INPUT THE CONFIGURATION 
7VGR-77 TfiRE ONLY 

I NP UT RAD I US » WE I GHTT » XBART » YBRRT > SGMRWT » SIGMRT » SGMCRL 
?£4. 1£7? 0. t 0. » 0. > 0. j 0. j . 0191 
INPUT X0*Y0» THETR IN BG. 

?-l U. 4333 0£s -6. 0£1080»-£9. 9734 
INPUT Ria»J> RDW BY ROW 

?37. 9QQ» 37. 934* 37. 979? 37 . 830 » 37 . 881 * 37. 76£j 37. 988* 37. 866* 37. 943 
INPUT 1 OR 2 FOR ONE OR TWO TRBLES OF L/C RERBINGS 
?1 


1. MERSURED C.M. 
CONFIGURATIONS 


TEST NO. 1 
RNB MASS 
VGR-77 TRRE ONLY 


INPUT 

MfiSS READINGS <KG!> 


900 

37.934 

37.979 




330 

37.881 

37. 76£ 




988 

37.866 

37.943 




RADIUS 

<CM> = 

£4. 1H70 




SGMCRL 

<KG> : 

. 0191 




OUTPUT 





XBRR 

<CM> : 

. 0048 

XBAR 

<IN) s 

.0019 

YBRR 

(CM) s 

0096 

YBRR 

<IN) s 

0038 

SIGMR 

<CM>: 

. 0035 

SIGMR 

WEIGHT 

<IN>: 

.0014 

NASS 

<KG) s 

113.6943 

<LB^ s 

£50.6530 

SGMRW 

<KG) : 

. 0£65 

SGMRW 

<LB) : 

.0584 




TORE 


INPUT 


XBRR 

<CM> s 

0.0000 

YBRR 

£CM> : 

0. 0000 

SIGMR 

<CM> : 

0. 0000 

MfiSS 

<K6> : 

0. 0000 

SGMRW 

<KG> i 

0. 0000 


3. MERSUREB VSR-77 TfiRE ONLY C.M. RNB MASS 
OUTPUT 


XBRR 

<CM> : 

. 0048 

KERR 

<IN> : 

. 0019 

YBflR 

<CM>: 

0096 

YBRR 

<IN> : 

0038 

SIGMR 

<CM> : 

. 0035 

SIGMR 

<IN> : 

. 0014 

MRSS 

<K6>: 

113.6943 

WEIGHT 

<LB> : 

250.6530 

SGMRW 

<K'G> s 

. 0265 

SGMRW 

s 

. 0584 


> INPUT TEST NO. 

72 

INPUT 1 IF UNITS RRE METRICS 2 IF UNITS ARE U.S. CUSTOMARY 

71 

INPUT THE CONFIGURATION 

7VGR-2 PARTIAL S/C DRY W/Q MIR IS 

INPUT RAD I US ? WE I GHTT ? XBART ? YBRRT ? SGMAWT ? S I GMAT ? SGMCRL 
724. 127 Oj 1 13. 6943? . 0048? 0096? . 0265? . 0 035? . 0650 
INPUT X0?Y0? THETA IN DG. 

7- 1 0 . 4338 02 ? -6 . 02 1080? -29 . 9784 
INPUT R1(I?J) ROW BY ROW 

7279.505? 277. 727? 244. 869? 276. 698? 243. 813? 27S. 577? 245. 354? 280. 008? 277.965 
INPUT 1 OR 2 FDR ONE DR TWO TABLES OF L/C READINGS 

72 

INPUT R£ ( I ? J) ROW BY ROW 

726 0 . £48 ? 253 . 286 ? 288 . 449 ? 252 . 334 ? £87 . 392 ? £59 .430? 288 . 956? 260 . 867 ? 253 . 609 


TEST NO. 


1. MEASURED C.M. RND MASS 

CONFIGURATION* VGR-2 PARTIAL S/C DRY W/O MIRIS 


INPUT 

MASS READINGS (KG) 



CASE 1. 



CASE 2. 


279.505 

277.727 

£44.869 

£60.248 

253.286 

288.449 

£76. 698 

£43.813 

278.577 

£52.334 

£87.392 

259.430 

245.354 

£80. 008 

£77.965 

288.956 

£60.867 

£53.609 


RADIUS (CM); 24.1270 

SGMCRL (KG) : . 0650 


OUTPUT 


XBRR 

(CM) s 

-. 0675 

XBRR 

(IN) : 

-. 0266 

YBAR 

(CM) s 

.9839 

YBRR 

(IN) : 

,3874 

SIGMA 

(CM) : 

. 0417 

SIGMA 

(IN)! 

. 0164 

MASS 

(KG) : 

801.5145 

WEIGHT 

(LB): 

1767. 0365 

SGMAW 

(KG) s 

. 0860 

SGMRW 

(LB)! 

- 1895 




OF 



x 




35 


£. TRRE 


INPUT 


XBRR 

<C.M> : 

.0048 

YBRR 

<CM> : 

0096 

SIGMR 

<CM> i 

.0035 

MRSS 

<KG> * 

113.6943 

SGMRW 

{KG!) ■ 

. 0265 


3. MERSUREB VGR-2 PRRTIRL S.'C DRY U)/0 MIRIS C.M. RND MRSS 
OUTPUT 


XBRR 

<CM> : 

-. 0786 

.YBRR 

<IN> ■* 

-.0310 

YBRR 

<CM!> * 

1.1466 

YBRR 

<IN> : 

.4514 

SIGMR 

<CM> s 

. 0485 

S I bMR 

<IN) : 

. 0191 

MRSS 

<KG> : 

687.8202 

WEIGHT 

<LB> : 

1516.3835 

SGMRW 

<K6> : 

.0900 

SGMRW 

<1_B> : 

. 1933 
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Appendix F. DERIVATION OF THE MEASUREMENT METHOD 


There is an advantage in using only one load cell successively to measure 
C.G. location. If a single load cell is employed, successively, in three 
known positions of a statically determinate fixture, with inert structural 
elements in the other two cells, the accuracy of the determination does not 
require an accurate knowledge of load cell calibration, but rather requires 
only that the load cell readout be linear with the applied load. An extension 
of this concept to a system incorporating three load cells and three succes- 
sive sets of measurements in which each load cell senses each of the three 
reactions gives an improved accuracy. 

Load cell calibrations are needed only to obtain a measurement of space- 
craft weight, hot for the determination of the center of gravity. 

1. Equilibrium 

Fig. F.l defines the configuration, the notation, and the sign conventions 
chosen for this analysis. 

At the outset,, it Is assumed that the manufacturing capabilities can 
provide a load cell placement accuracy that removes planform geometry from 
consideration as a significant contributor to the error. Thus the principal 
concern relates to errors associated with load cell readouts. Assume that the 
reactions due to the weight of the spacecraft are R^, Rg and R^. The weight 
and C.G. platform axes are designated as X* and Y'. For the present, it is 
assumed that the spacecraft Z axis coincides with the platform Z’ axis through 
the centroid of the equilateral triangle defined by the reaction points. 
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Fig. F.l. System Geometry and Notation 


Thus, in the spacecraft system, 

X - X f - | : Y = Y* 
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The total weight is 


w - ^ + r 2 + r 3 


From the moment equilibrium about the Y' axis 


=E i(i r ) 


_ 3 "-i 

yl - 2 ■*» 

X 2 W 


From the moment equilibrium about the X 1 axis. 


Hf ' ‘I ( R 2 - R 3 ) 1 


2. Error analysis based on a single set of measurements. 

In Eqs. (1) to (3), the reactions R,. are taken to be actual values, 
without error. Let the load cell readouts be defined as 


P. = k.R. (4 

3- i J 

where k^ is a nondimensional constant. Then using P^ in place of 

in Eqs. (1) and (2), and expressing errors in differential notation leads to 


— ■ ' — ' HP — 1 

3x_dr 3x l + w 3x_ £W 

r 3r r 1 8W W 


where 


P ™_!!i + P i!L ^2. + P 

i 3p i p i 2 Sp 2 p 2 3 ap 3 p 3 




dr - 1 dP l ( P 1 + Zl ^jL + ll 

dx = x — + x ~y~ - x ^ + w p 2 ^ W p 3 


As finite differences, considering that AP. may be either positive or negative. 


the "worst case" gives 


, Ar , /, P l\^l + !2.!!2 !3^3 

CAx >MAX ‘ X ~ + l 1 " J P. » P, W P 3 


Similarly, it can be shown that 


(Ay' ). 


, ? -(a, + ^ + ^ + ^) 

\ r p x p 2 P 3/ 


+ 2 1 1 W P 2 W P 3 


If it is assumed that 


p. = P; AP, « AP = AP 
3 1 J 

r =■’ _ r\/5" 


then Eqs. (5) and (6) become 


__t i 2 AP 

(Ax ^MAX ~ 2 Ar + IT r 


a/3 AP 


CAY ■’MAX “ 3 P r 




3. Examination of AP/P, 


The K in Eq, (4) are nondimens ional constants, and can be written 

k. = E. S. V. 
x x i x 


( 9 ) 


where V. is the voltage readout from the load cell circuit, S. is the basic 
x x 

sensitivity of the load cell in lbs/volt per volt, and E^ is the voltage frc 
the bridge power supply. Dropping subscripts, we have, in general 


P = ESVR 


dP _ E 3P dE S 3 P dS V 3 P dV . R 3 P dR 

P P 3 E E P3S S + P 3 V V P3R R 

, * dR = o 

Since, by definition, R is the actual load . . ^ 


( 10 ) 


ap _ ae As av 

P E S V 


( 11 ) 


Regulated D.C. power supplies are available with a voltage stability 
AE/E of 0.005% over the test period. If a digital voltmeter is used for read- 
out, there is an error (AV/V)^ of plus-or-minus one digit in the last place 
during (1) bridge balancing, and (2) during load readout. Considering that 
these errors can occur during load cell calibration and again during weight- 
and-C.G. measurement, the worst-case error where the calibration is relied upon 


when readout to the nearest tenth-pound is provided. The other part of error 
(AV/V)^ relates to the reproducibility of the voltmeter from one measurement to 
another s rather than to the absolute calibration in volts. For the purposes of this 
analysis, the reproducibility is assumed to be 0.02% as an upper limit. Thus 

* 

= 0.00005 + ^ + 0.0002 + 

Assume P = 500 lbs, then 

= 0.00105 + jr- (12) 

The basic sensitivity, S , can be deduced from the relation 


where 


S 


02 

EVW 


(13) 


£2 = resistance in ohms of a shunt calibration resistor 
across one leg of the load cell bridge 

C = "equivalence factor," lbs/ohms 

W = calibration weight, lbs 


To clarify Eq. (13) it is noted that a standard calibration procedure 
is as follows: 

(1) balance the bridge with no load 

(2) load the cell with a calibration weight, W , near the rated 
capacity of the cell 

(3) adjust the sensitivity of the readout box until the digital 
display coincides with the weight, W . 
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(4) unload, and check zero. 

(5) repeat these steps until zero load and calibration load 
readouts agree with actual loading 

(6) insert shunt resistor and note "equivalent load" 

0 (7) proceed with calibration in load increments to obtain a 

check of linearity. 

If the readout, P , is nonlinear with calibration weight, W , the 
calibration curve may be used in lieu of a constant sensitivity, S . All 
that is required for accurate measurement is that the calibration curve is 
reproducible . 

An error analysis of S , based on Eq. (13) is not presented here, as 
it will be shown that it is of no practical concern in C.G. determination, 
‘"although it is for the weight determination. 

4. A rudimentary C.G . measurement scheme . 

Asstime that the test fixture uses only one load cell, with inert 
elements at the other two apices of the equilateral triangle. Assume that 
three separate measurements are made, with elements of the fixture rotated 
so that the load cell is subjected, successively, to the reactions R^, R 2 
and Furthermore, assume that the load cell sensitivity, S , is not 

independent of load, as depicted in the sketch below. 
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Then, with 


and Eq. (12) becomes 


P l“ P 2” P 3' i f ” 0 


AP 

^ w 0.00105 

Experience with load cell calibrations shows that nonlinearities are very 

small, and that the assumptions here are valid within measurement demands. 

Let us now assume that r = 24 in., Ar = .004 in. Then, from Eqs. (7) and (8) 
(AKJjfc x " 0*021 in. 


(AY ) max = 0.015 in. 


/(AX’)^ + (Ay')^ = 0.026 in. 


This is a very good accuracy. 


5. The implemented system. 

The system implemented will use three load cells, ' instead of one. Here 

P ij = k i P j (14) 

where the subscript, i , refers to a particular load cell and the subscript, j , 
refers to a particular reaction. By defining the average G.G. coordinates as 


I Cx’j + x 2 + x' 3 ) 


I (y'i + i’ 2 + y’ 3 ) 


a differential analysis reveals that the maximum errors are the same as those 


given by Eqs. (7) and (8) 
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Oth er options . 

The sections above discussed only one option in the assessment of errors 


in C.G. determination. This option entailed the use of instrumentation giving 

load cell readout to the nearest tenth-pound ani a test procedure of three 

sets of measurements such that each of the three load cells measures each of 

the three spacecraft reactions. Three other options are presented here for an 

absolute calibration AS _ qq 2 

S 

Case 1. One pound readout? 3 load cells with one set of measurements 

Then Eq. (12) of the reference becomes 

AP 4 

— = .00005 + .002 + .0002 + g-g 

= .00725 

From Eqs. (7) and (8) 

(“'W =-118 in. 

(“'Jm ’ - 100 

4r ilAX = in. 

Case 2. Five- tenth pound readout; 3 load cells with 1 set of measurement. 


— = .00005 + .002 + .0002 + .0025 
P 

= .00475 
) mAy - *0^ i n * 

U*'W= - 0658 to - 

( ir'L„ - .102 in. 


Case 3. One- tenth pound readout; 3 load cells with one set of measurements 


— = .00005 + .002 + .0002 + .0005 
P 

= .00275 

( a *')max = * 046 ±n * 

( A *’W = * 038 in - 

( Ar * -MAX = *° 6 ° in - 
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Appendix G. STATISTICAL BEST FIT 


SUMMARY 

Weight and center of gravity measurement technique has been established 
as discussed in Sec. 2 of this report. Six locations of the apparent C.G. 
of the spacecraft are determined by making three 120° load cell rotations and 
one 180° spacecraft rotations to eliminate systematic errors. The present 
analysis finds two best fit equilateral triangles to the six points in a least 
square sense in order to improve the C.G, accuracy and give an error estimate. 

The centroid of the two equilateral triangles is taken as a ’’better" 

C.G. and the standard deviation is calculated. 


ANALYSIS 

This analysis is aimed at minimizing instrumentation errors. Sec. 2 of 
this report shows that for each spacecraft rotation the apparent C.G. locations 
for three 120 load cell rotations should form an equilateral triangle in 
absence of readout errors. The present analysis shows that the mean value is 
actually a best fit. The standard derivation o for the obtainable accuracy 
is determined by this analysis and indicates what percent of the errors can be 
attributed to instrumentation by statistical estimate. 

Six data points (x^y^ , i = 1 6 are given by Eqs. (1), (2), (3) of 

Sec. 2. By rotating a vector whose x,y components are a and b through 
two 120° angles, the vertices of an equilateral triangle are obtained, located 
close to one set of three data points corresponding to one spacecraft position. 
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•if. 





If (x,y) Is the new C.G. and y,e are the vector components from this C.G. 

to the centroid of either triangle (see Fig. G.l) the coordinates of the 

vertices of one triangle can be written as (x i y i a. , y - e + b^) where 

i = 1 9 . 3 and (x - y + y-fe-f b^) where i = 1, ... 3 for the other, 

2 

Let S be the mean square error. 


2 r i i j 

Zj x °i - (x + Y + a ± ) + ^ y°. - (y - e + b ) 

i=ll J i=l L 

3, r -|2 3, P 

+ A. *“ 1+3 - - -r + a ) + 5 7° 1+3 - (y + s 

1=1 L J i=l I 


1 * Jr 

* 2 A 

J 1=1 L 


^3 (y + £ + b_^) ( 1 ) 


2 2 — 2 — 

The mean square S is minimized by equating the partials 3S /3x, 3S /3 y 
2 2 2 2 

3S /3y, 3S /3e, 3S /3a and 3 S /3b to zero. This set of equations is solved to 


o n 

1 V o - 1 V o 

" e A x i’ y = ?A y i 


'3 6 1 3 6 

1 y o y o 1 y a y □ 

"6 Z x i -Z x i > e = 6 Zy^Zvi 

D 1=1 1 1=4 i=l 1 1=4 


1 „ . o , o . \ o r~ , o o,o o. 

" 12 3Cx 1 + x 4 } ~ A x i + ^ (y 2 " y 3 + y 5 ~ y 6 } 


b "Ij[ 3( A + - 5, y\-^ 

1=1 


/ o O t O O v 

(x 2“ x 3 + x 5“‘ x g) 


It is seen that (x,y) is the mean of the data points. The standard devia- 
tion a due to errors made in the measurements can now be estimated from Eq. (1) ♦ 
Noting that there are n =6 points, an unbiased estimate of the 
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( 6 ) 


standard deviation is obtained by dividing 


by n -1=5. 



(vector C.G.) 


Standard Deviation for Weight 

The values W^, W^ and refer to the weight of spacecraft and tare 
for one spacecraft position and W^, W^ and Wg for a 180° rotation of the 
spacecraft. Recalling Eqs. (4) and (7) in Sec. 2, the standard deviation cr^ 
for the total weight can be calculated from 


a w = 


~ W) 2 + ... + CW 6 - W) 2 2 

. 5 + ° CAL 


4/2 


(7) 


where W is the mean value of the U4s and is the standard deviation 

from the calibration curve of the load cells. 


Standard Deviation for Spacecraft 

Given tare data W T> Y^, and o Wt we calculate \/ c ^ s / c and w s / c 
from Eqs. (8), (9) and (10) in Sec. 2, and also the standard deviations for: 

(a) Spacecraft C.G. 


a 


_ w 

s/c " w s/c 


a 


( 8 ) 


b) Spacecraft Weight 



0EIGWALPM5BK 

of poc® $ UALU 


( 9 ) 


A sample output and the program listing are in the Appendix E and A, 


respectively. 
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Appendix H. THE WEIGHT AND THE CENTER OF GRAVITY ERROR ESTIMATES 


1. SUMMARY 

It is noted that the effect of the spacecraft C.G. offset (from the 
center line of the load cell assembly) on the spacecraft C.G., error estimate 
is eliminated by taking two sets of measurements with the spacecraft rotated 
180° between measurements. The obtainable accuracy on the spacecraft C.G. 
from the Voyager spacecraft weight and C.G. measurement fixture, with more 
realistic weight estimates and a more conservative tare C.G. location error 
estimate than Sec. 2, is estimated to be ±0.017 in. And, the corresponding 
accuracy in the spacecraft weight is estimated to be ±4.5 lbs. 

2. ANALYSIS 

From Eas. (5) and (6) of Appendix F. 


- -’Ur +n V P 1 + P2 ^ + !l ^3 

(Ax >MAX = * [— + C1 " l? "W P 2 W P 3 


(Ay ), 


. AP. AP_ AP, 

Al + _1 +_! + _ 1 

r ^ P P P 

r 1 2 3 


12^2 A 

W P 2 W Pj 


where P^ + Pg + P^ = W 


These equations are simplified by assuming that 


P i * F j = p 5 x « f and y' « 


Furthermore, AP. a AP. = AP, since the three load cells must have the 
i 1 

same bound. Therefore, 


, A -\ _ Ar , 2 AP 

^ Ax \l&X ~ 2 + 3 P r 


_ \/3~AP 

^ Ay ^MAX 3 P r 


The assumptions above can be relaxed with no change in Eqs, (7) and (8) 
other than 3P + W. The two sets of measurements taken with the spacecraft 
(with C.G. offset of x,y) rotated 180° (not part of the derivation of Appen- 
dix F.) therefore eliminate the effect of C.G. offset of the spacecraft on the 
error analysis. 


„ » 


r . — 


For, Position 1 x = + x 


(180° from 1) 




r — 


Position 2 x = — - x 


— * = -T 


From Eq* (5) above 


t 


Ar ( Wx 1 \ AP 1 — ' ^2 AP 3 

i £ i t + - x J“ir + !£ iir + x i"r 


and from Eq. .(2) of Appendix F. 


_ 3 1 

X 1 2 W r 


Therefore, 


Ax 'i s 






iP 2 , 4P 3 
-r + x i — 


or 



Ar 

1 



AP 

W 


Similarly, 


Ax’ 2 £ 




3 

2 r + X 



AP 

W 
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Then 


, A , + iP 

2 2 r + \2 r + 2 7 ~W 


Therefore 


or 


fix' S^+ 2r^| 


_T 2 

Ax ^ g ^ “p r which is the same as Eq. (7) with W « 3P 


From Eq. (6) above 


or 


Similarly 


(^ + f) + # (2 f ) 




a?' s 7' (^s. + m.) + r Jl to 

y 2 y 2 > r »/ V3 W 


Then 


t _» t t 

_> Ay + Ay y . + y 

A v = — S — i 

7 2 2 V r ' W 




Ap 

w 


Therefore 


A7' s 


or 


— V3 AP 

Ay ^ — j — r which is the same as Eq. (8) with W w 3P 
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3. ERROR ESTIMATE 


a) C.G. location of spacecraft with Tare, and C.G. location of spacecraft. 
From Eqs. (11) and (12) where 3P is replaced by W of Sec. 2 

AX $ —• + 2r — + Ax 
2 wo 

AY ( r /3 ^ + iy 
W o 

Where Ax Ay are concentricity error of the spacecraft longitudinal 
o’ J o 

axis relative to the vertical reference axis of the weight and C.G. fixture 
as mentioned in Appendix F. 

Let 

Wg /C = 1650 lbs , W T = 300 lbs 

r = 9.50 in. Ar = 0.002 in. 

Ax = Ay = 0.0035 in. 
o J o 

and AP = AP X + AP 2 + AP 3 

The readout error AP is a random error consisting of three parts, 

AP^, AP^ and AP„ where 

i. AP^ is that part of the error due to the readout drifting during 

calibration. Its value has been determined to be 0.10 lb from the 
load cell calibration curves over the range of 0 ■> 960 lbs. for all 
three load cells. Only one reading among a total of 117 gives a 
AP^.>- 0.10 lbs. and was discarded as an anomaly. 
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ii. AP£ Is that part of the error due to the readout resolution. Its 
value has been chosen to be 0.1 lbs. 

Ui . AP^ is that part of the error due to the readout fluctuation on 

account of the test site environment which includes ground vibration, wind 
draft, etc. Its value was not known so a range from 0 + 0.50 lbs 
was assumed. 

Then the C.G. error was estimated giving the following table: 


;£.-3 




•u 1 


1 » 

«• 


a 


AP 3 

(lb) 

AP 

(lb) 

AX 

(in.) 

AY 

(in.) 

ACG S/C + T 
(in.) 

^S/C 

(in.) 

A *s/c* 

(in.) 

4CGs /c l 

(in.) 

0 

0.20 

0.0065 

0.0052 

0.008 

0.0077 

0.0062 

0.010 

0.10 

0.30 

0.0074 

0.0060 

0.010 

0.0088 

0.0071 

0.011 

0.20 

0.40 

0.0084 

0.0069 

0.011 

0.0099 

0.0082 

0.013 

0.30 

0.50 

0.0094 

0.0077 

0.012 

0.0111 

0.0091 

0.014 

0.50 

0.70 

0.0113 

0.0094 

0.015 

1 

0.0134 

1 

0.0111 

0.018 


*Calculated from Eqs. (13) and (14) of Sec. 2. 

b) The error of weight measurement has been estimated by the following 
equation 

AW AX , AP 

~w = T + “? 


or 

where W M 3P 

AX 


AW = ^ W + 3AP 
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and 


is the load cell calibration accuracy, assumed, to be 0.1%. 
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